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Abstract 

, We present the predictions of a model for proton-proton total cross- 

,_i | section at LHC. It takes into account both hard partonic processes and 

soft gluon emission effects to describe the proper high energy behavior 
^ . and to respect the Froissart bound. 

^ <— 1 1 Introduction 

A reliable prediction of the total proton-proton cross section is fundamental to know which will 
be the underlying activity at the LHC and for new discoveries in physics from the LHC data. In 
this article, we shall describe a model [1] [2] for the hadronic total cross section based on QCD 
minijet formalism. The model includes a resummation of soft gluon radiation which is necessary 
, to tame the fast high-energy rise typical of a purely perturbative minijet model. It is called the BN 

model from the Bloch and Nordsiek discussion of the infrared catastrophe in QED. In the first 
section, results are presented concerning the behavior of the QCD minijet cross section. It will 
then be explained how this term is included into an eikonal formalism where infrared soft gluon 
emission effects are added. The last section is devoted to the link between the total cross-section 
r*** ' asymptotic high energy behavior predicted by our model and the model parameters. This relation 

also shows that our prediction is in agreement with the limit imposed by the Froissart bound. 

o : 

J> . 2 Mini-jet cross section 

Hard processes involving high-energy partonic collisions drive the rise of the total cross section 



[3]. These jet-producing collisions are typical perturbative processes and we can describe them 
through the usual QCD expression: 



Ofct(. s iPtmin) = dp t dxi dx 2 x fi\A(xi,p 2 )f j \ B (x2,Pt)—T- 

J P train 

Jipj/s J4 P 2 t /(x lS ) .... ap, 

(1) 

with A, B = p,p. This expression depends on the parameter ptmin which represents the min- 
imum transverse momentum of the scattered partons for which one allows a perturbative QCD 
treatment. Its value is usually around rj 1 — 2 GeV and it distinguishes hard processes (that 
are processes for which a perturbative approach is used) from the soft ones that dominate at low 

' speaker 



energy, typically for i/s < 10 4- 20 GeV, i.e, well before the cross-section starts rising. The 
Minijet expression also depends on the DGLAP evoluted Partonic Densities Functions f^ A for 
which there exist in the literature different LO parameterizations(GRV, MRST, CTEQ [4]). We 
obtain an asymptotic growth of <x, e t with energy as a power of s. As shown in figure [T] the value 
of the exponent depends on the PDF used and one has 

„GRV ^ o 0.4 MRST ^ n.3 CTEQ ^ n.3 
a jet ~ s a jet ~ s a jet ~ S ■ 

This result can be derived by considering the relevant contribution to the integral in (Q} in the 
>> ptmin limit. In this limit, the major contribution comes from the small fractions of 
momentum carried by the colliding gluons with xxp << 1. In this limit we know that the relevant 
PDF's behave approximately like powers of the momentum fraction x~ J with J ~ 1.3 [5]. 

From the previous consideration and noting that — ^ — oc ^ we obtain from £T|) the following 
asymptotic high-energy expression for aj et : 

1 
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(2) 



t min_ 

The dominant term is just a power of s and the estimate obtained for the exponent e = J— 1 ~ 0.3 
is in agreement with our previous results. We now need to understand how to incorporate into 
a model for the total cross section this very fast rise at very high energy, which is present in 
the perturbative regime. Firstly it is important to note that <jj e t is an inclusive cross section and 
therefore contains in itself a multiplicity factor, linked to the average number < n > of partonic 
collisions that take place during the hadronic scattering. We can approximate the energy driving 
term at high energy [6] < n > as 

< n >ss a j e t ■ A, (3) 

where A is a function representing the overlap between the two hadrons. 

Now we can derive an expression for the total cross section as a function of < n >. Assum- 
ing that the number of partonic collisions follows a Poisson distribution, since each interaction is 
independent from the other, the probability of having k partonic collisions is: 

< n p~ < n > 

P(k,<n>)= fc , • (4) 

The average number of partonic collisions should depend on the energy and on the impact pa- 
rameter b relative to the hadronic process < n >=< n(b, s) >. From the previous expression it 
is possible to obtain the inelastic hadronic cross section: 



^inelastic = [ dH £ P(k, < n(b, »)>)=/ d% [l - e~ <n ^ S > 
J fc=l J 



(5) 



which is the usual eikonal expression if we consider the link between < n(b, s) > and the eikonal 

X{b,s): 

< n(b, s) >= 2Imx(fe, s). (6) 
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Fig. 1: minijet cross section for different input parton densities. 



3 Eikonal model 

The eikonal representation allows to implement multiple parton scattering and to restore a finite 
size of the interaction. Neglecting the real part of the eikonal function, an acceptable approxima- 
tion in the high energy limit, the expression for the total cross section is 

1 - e - nib ' s)/2 ] . (7) 

The average number of partonic collisions receives contributions both from hard and soft physics 
processes and we write it in the form 

n(b, s) = n soft (b, s) + n hard (b, s), (8) 

where the soft term parameterizes the contribution of all the processes for which the partons 
scatter with p t < ptmin- It is the only relevant term at low-energy and it establishes the overall 
normalization, while the hard term is responsible for the high-energy rise. From (O, we approx- 
imate this term with 

nhard(b, s) = A(b, s)a jet (s), (9) 

where the minijet cross section drives the rise due to the increase of the number of partonic 
collisions with the energy and A(b, s) is the overlap function which depends on the (energy 



a t ot = 2 d 2 b 



dependent) spatial distribution of partons inside the colliding hadrons. In some older models 
[6] a simpler factorized expression for n(b, s) was used, with the overlap function depending 
only on b. However, when up-to-date realistic parton densities are used, such impact parameter 
distributions, inspired by constant hadronic form factors, led to an excessive rise of a to t with the 
energy. In our BN model we include an s-dependence in the overlap function that has to tame 
the strong growth due to the fast asymptotic rise of Oj e t [2]. 

We identify soft gluon emissions from the colliding partons as the physical effect respon- 
sible for the attenuation of the rise of the total cross section. These emissions influence matter 
distribution inside of the hadrons, hence changing the overlap function. They break collinearity 
between the colliding partons, diminishing the efficiency of the scattering process. The number 
of soft emissions increases with the energy and this makes their contribution important, also at 
very high energy. The calculation of this effect uses a semiclassical approach based on a Block- 
Nordsieck inspired formalism [7], the basic assumption of this technique is that all emissions 
are independent from each other, so the number of gluons emitted follows a Poisson distribu- 
tion. Thereof one obtains a distribution of the colliding partons as function of the transverse 
momentum of the soft gluons emitted in the collision, i.e. 



d 2 P(K ± ) = d 2 K ± j^ J d 2 b e iK ^ h - h ^ , (10) 
the factor h(b) is given by 



iif>; />/,,•;/.•); i ( '''' '•; j'^ X) lf ,m (%nW|[i- '" 

" m,n=colors 

where d 3 n g (k) is the distribution for single gluon emission in a scattering process and it is linked 
to the QCD current responsible for emission. 

We have proposed to obtain the overlap function as the Fourier transform of the previous 
expression of the soft gluon transverse momentum resummed distribution, namely to put 

/r/ 2 P(\C i p -h(b,q ma x) 

with 

,16 f Qmax aJk?) dk t , 2o maxr , 
h(b,q max ) = — - j}l±L -t "log^f^ 1- Jo (k t b )], (12) 
3 Jo n kt kt 

this integral is performed up to a maximum value which is linked to the maximum transverse 

momentum allowed by the kinematics for a single gluon emitted, q ma x [8]. In principle, this 

parameter and the overlap function should be calculated for each partonic sub-process, but in 

the partial factorization of Eq.© we use an average value of q ma x obtained considering all the 

sub-processes that can happen for a given energy of the main hadronic process [2] : 

" max(s) = v 2 yJWJWLJ^^) ' (13) 



with z m i n = Apf min / \sx\X2). Notice that consistency of the calculation requires that the PDF's 
used in Eq. (PT3T ) be the same as those used in Oj e t- In Figj2]are presented our results for q max as 
function of y/s using ptmin = 1.15 GeV. 




Fig. 2: Results for q ma x using different input partem densities with ptmin =1.15 GeV. 



The integral in (fT2l has another relevant feature, it extends down to zero momentum val- 
ues, and to calculate it we have to take an expression of a s different from the perturbative QCD 
expression which is singular and not integrable in (fT2l ). We use a phenomenological expres- 
sion [9], which coincides with the usual QCD limit for large k t , and is singular but integrable for 

2 127T p 

as[kt) = 33-2iV / ln[l+p(^)2p]- (14) 

This expression for a s is inspired by the Richardson expression for a linear confining potential 
[10], and we find for the parameter p that 

• p < 1 to have a convergent integral (unlike the case of the Richardson potential where 
P = 1) 

• p > 1/2 for the correct analyticity in the momentum transfer variable. 

Fig]3][l] shows our predictions, obtained for the total cross-section using a set of phenomeno- 
logical values for ptmin and p, and varying the parton densities. We also make a comparison with 
data and other current models. 




Fig. 3: Results from our total cross-section model (for different parton densities) compared with data [11] and with 
other models [12], 



4 Restoration of Froissart Bound 

The Froissart Martin Bound [13] states that otot cannot rise faster than a function which is pro- 
portional to log 2 (s). In order to see that in our model this bound is respected, we approximate 
our total cross section at very large energies as 



a tot w 2tt I db 2 



1-e 



-n hard (b,s)/2 



(15) 



with rihard(b, s) (Tj e t(s)Ahard(b, s). We then take for aj e t the asymptotic high energy expres- 
sion: 

with o"i =constant and e ~ 0.3 — 0.4. Being Ah ar< i(b, s) oc e~ h<yh,s \ we can consider in (fT2l ) the 
infrared limit k t —> where the integral receives the dominant contribution. In this limit we have 

2^ _ / IV \ 



apart from logarithmic terms. Then, with h(b, s) oc (bA) 2p [2] (again apart from logarithmic 
terms), we have 

A ha rd(b) oc e 



-(bA) 2 P 



and from this expression 

n hard = 2C(s)e-^ 2p , 
with C(s) = ^9p- ( G *y± ) ■ The very high energy limit of Eq. (fT5l) then gives 

f^[l-e-^- (6A)2P ] 
J o 

The asymptotic growth of atot in our model depends on the parameter e which fixes the asymp- 
totic rise of the minijet cross section, and on p which modulates the infrared behavior of a s . 
Notice that 1/2 < p < 1 and thus this approximated result links the restoration of the Froissart 
bound in our model with the infrared behavior of a s . We can now understand why a knowledge 
of the confining phase of the strong interaction is necessary if we want to restore the finite size 
of the hadronic interaction. 
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